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Eccentricity fluctuations make flow measurable in high multiplicity p-p collisions 
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Elliptic flow is a hallmark of collectivity in hadronic collisions. Its measurement relies on analysis 
techniques which require high event multiplicity and could be applied so far to heavy ion collisions 
only. Here, we delineate the conditions under which elliptic flow becomes measurable in the samples 
of high-multiplicity (dN c h/dy > 50) p-p collisions, which will soon be collected at the LHC. We 
observe that fluctuations in the p-p interaction region can result in a sizable spatial eccentricity 
even for the most central p-p collisions. Under relatively mild assumptions on the nature of such 
fluctuations and on the eccentricity scaling of elliptic flow, we find that the resulting elliptic flow 
signal in high-multiplicity p-p collisions at the LHC becomes measurable with standard techniques. 
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In high-energy hadronic collisions, all projectiles, be it 
protons or nuclei, have finite spatial extension and thus 
collide generically at finite impact parameter. This fi- 
nite impact parameter, as well as event-by-event fluctu- 
ations, both result in an azimuthally asymmetric shape 
of the interaction region in the plane transverse to the 
beam direction. If particle production is not governed 
solely by independent local processes, it can be modified 
collectively by the azimuthally asymmetric spatial gradi- 
ents present in the collision geometry. In this case, the 
spatial eccentricity e of the interaction region results in 
a momentum anisotropy of single inclusive particle mo- 
mentum distributions dN in the azimuthal angle 0, which 
can be characterized by the harmonic flow coefficients v n 
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In general, the v' n s depend on transverse momentum 
and rapidity y. The measurement of finite harmonic flow 
coefficients v n is widely regarded as one of the most direct 
dynamical manifestations of collective behavior in high- 
energy hadronic collisions 0, Q . 

The measurement of flow coefficients v n is complicated 
by the fact that it requires information about the az- 
imuthal orientation of the reaction plane. As we 
recall below, whether this reaction plane is experimen- 
tally accessible depends on both, the strength of the flow 
signals v n and the multiplicity n mu it in the phase space 
window in which v n is determined. This is so since mi- 
croscopic dynamics such as dijet production is a source 
of azimuthal asymmetry uncorrelated to 4>r and needs 
to be disentangled from the collectively preferred macro- 
scopic motion of all particles, characterized by v n . Up 
until now, a sufficiently large signal strength at suffi- 
ciently high event multiplicity has been observed only 
in heavy ion collisions, and measurements of flow coeffi- 
cients v n have not been reported for p-p collisions. It is 
the main purpose of this Letter to provide generic argu- 
ments for why at the LHC collective flow coefficients v n 
may become experimentally accessible for the first time 
in proton-proton collisions and to discuss how this mea- 



surement would impact the understanding of soft physics 
both in p-p and in heavy ion collisions. 

There are several, mutually compatible standard tech- 
niques to discriminate collective flow coefficients from 
non-flow effects 0-0] • Here, we focus on the cumu- 
lant analysis [3] which is based on the observation that 
azimuthal correlations between arbitrary particle pairs 
in an event, 
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do not arise solely from global correlations with (j>R 
(which give rise to v n ), but that they can also have 
other, microscopic origins 5. Since the non-flow effects 
parametrized by S are, by definition, correlated at most 
to a small subset of all particles in the event, they are 
parametrically suppressed by one power of n mu i t . Hence, 
the elliptic flow coefficient v-i can be reliably extracted 
from two-particle correlations if «2{2} > 1 / ^/n mu it . The 
assumption 3> S that measured azimuthal two-particle 
correlations are dominated by a global correlation with 
(j>n can then be tested systematically by studying higher- 
order particle correlations. For instance, one can analyze 
four-particle correlations (4th order cumulants) which re- 
quire only a signal size ^{4} > 1/ (n mu i t ) 3 ^ 4 . By going 
to higher cumulants, one achieves at best a sensitivity 
v-2 > 1/^muit- In general, the cumulant analysis allows 
one to disentangle the collective flows v n from non-flow 
corrections by systematically exploiting the different mul- 
tiplicity scaling of flow and non-flow contributions in dif- 
ferent n p -particle correlations. 

With such standard techniques, the lowest order col- 
lective flow coefficients, in particular v±, vi and v^, have 
been characterized in heavy ion collisions unambiguously 
and over a wide range of impact parameter at center of 
mass energies ^/snn 

between 2 and 200 GeV @-[!]. In 
the following, we focus mainly on elliptic flow v-i at mid- 
rapidity. Odd harmonic coefficients vanish in this case by 
symmetry and i>2 is known to parametrize the dominant 
momentum space asymmetry. The size of V2 measured in 
experiments at the Relativistic Heavy Ion Collider RHIC 
are very large. For instance, for transverse momenta in 
the range px — 2 GeV, semi-peripheral Au-Au collisions 
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at RHIC result in more than twice as many hadrons pro- 
duced in the direction of the reaction plane, than orthog- 
onal to it (v 2 (pt ~ 2GeV) > 0.2). The transverse mo- 
mentum integrated v 2 reaches i> 2 — O(0. 1). Remarkably, 
the experimental praxis of -^-measurements indicates 
that the parametric bounds on v 2 {2}, v 2 {4}, ...w 2 {n} 
given above provide realistic numerical estimates for the 
feasibility of v-i measurements if one uses for n mu it val- 
ues of order of the charge multiplicity per unit rapidity 
dN c h/dy. We summarize this information in Table U for 
"mult = 30, 50 and 80. These values of n mu it are smaller 
than dNch/dy in sufficiently central heavy (Au) collisions, 
but they are comparable to the values in semi-peripheral 
collisions of lighter (Cu) nuclei. 
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> 0.18 


> 0.09 




tumult 


= 50) 


> 0.14 


> 0.05 


t>2 


(^mult 


= 80) 


> 0.11 


> 0.04 



TABLE I: Estimates of the minimal signal strength V2{n p }, 
which can be discriminated from non-flow effects in an n p - 
cumulant analysis based on n mu i t particles. 

The interpretation of elliptic flow measurements in 
heavy ion collisions relies on the observation that u 2 is 
correlated with the initial spatial eccentricity e of the 
transverse overlap region of the two projectiles [l| , 
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Here, averages are performed with respect to the mat- 
ter distribution right after the collision, and x 1 and y' 
denote the lengths along the main axis of an ellipsoid de- 
scribing this distribution. For sufficiently central heavy 
(Au) and lighter (Cu) ion collisions, it is found that 
«2 a e [l(J EH ■ Remarkably, this is a generic expecta- 
tion of fluid dynamic simulations of heavy ion collisions. 
Models of ideal dissipation-free hydrodynamics, which by 
construction describe collision scenarios of maximal col- 
lective flow, can account quantitatively for the size of the 
elliptic flow measured at RHIC [12h14| | and the expected 
dissipative corrections are anomalously small 15H18|. In 
conjuction with this interpretation, the observation of 
very large v 2 -signals is arguably one of the most far 
reaching discoveries of the RHIC heavy ion program. 

We now turn to the question whether elliptic flow may 
be measurable in p-p collisions at the LHC. In p-p col- 
lisions studied so far, one may explain the apparent ab- 
sence of an elliptic flow signal by pointing to the fact 
that the dN^/dy in these collisions is too low to make v 2 
measurable (see Table HJ. However, while Monte Carlo 
simulations for minimum bias dN^/dy distributions in 
,/s = 14 TeV p-p collisions peak at low values < 10 for 
the non-diffractive contribution, they show a pronounced 
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FIG. 1: Schematic view of region of hadron production may 
be located in the transverse overlap region of a central proton- 
proton and central nucleus-nucleus collision respectively. De- 
pending on the number and size of hadronically active regions, 
large eccentricities can result even in central collisions. 

high-multiplicity tail, typically reaching values as high as 
dN c h/dy ~ 60. Despite their model dependence, these 
simulations strongly indicate that abundant samples of 
high- multiplicity p-p events with dN^/dy > 50 will be 
measured at the LHC. Such a multiplicity is comparable 
to that reached in semi-peripheral (centrality class 40 - 
60 %) Cu-Cu collisions at JsWn = 62.4 GeV at RHIC, 
and for these latter collisions elliptic flow has been mea- 
sured. Whether elliptic flow is also measurable in high- 
multiplicity p-p event sample at the LHC then depends 
on the signal strength v 2 and on the relative strength of 
non-flow corrections in p-p collisions. 

To estimate the strength of the elliptic flow signal 
i> 2 , we now discuss the initial spatial eccentricity e of 
hadronic collisions. In the collisions of heavy (Au or Pb) 
ions, this eccentricity is determined solely by the trans- 
verse spatial overlap. More precisely, in a nucleus-nucleus 
collisions dN c h_/dy scales approximately with the average 
number iV part of participant nucleons, which scales with 
the area of the nuclear overlap. Therefore, selecting a 
multiplicity class in A-A amounts to selecting on impact 
parameter and determines the shape of the nuclear over- 
lap region. For large N paiIt , it is reasonable to make the 
smoothness assumption that the interactions between the 
^Vpart nucleons result in a homogeneous density distribu- 
tion within the area of the nuclear overlap (for illustra- 
tion, see right hand side of Fig. [T]) . If this assumption 
would carry over to p-p collisions, then the highest mul- 
tiplicity p-p collisions would be the most central ones, 
their spatial eccentricity would be close to zero, and so 
would be the flow signal u 2 oc e. Previous estimates of 
the magnitude of vi were based on this smoothness as- 
sumption fl9l [20l ] or on other methods [2lj and reported 
small, non-measurable values. 

However, sizable deviations from the smoothness as- 
sumption have been found in modeling lighter (Cu) ion 
collisions 1(J U|. In these systems, the relatively small 
number of nucleon-nucleon interactions results in event- 
by-event fluctuations of the density distribution which 
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can increase the initial spatial eccentricity significantly 
above its geometric estimate. The predicted eccentric- 
ity scaling of v 2 is only confirmed in Cu-Cu collisions, if 
these fluctuations are taken into account 1(| 

For p-p collisions, the A part -scaling of event multiplic- 
ity used in A-A does not apply. It remains a priori un- 
clear to what extent a cut on event multiplicity amounts 
to a cut on impact parameter. MC simulations of the 
underlying event in p-p do assume indeed a correlation 
between impact parameter and the number of indepen- 
dent partonic interactions Ampi, which determines event 
multiplicity and, depending on model assumptions, may 
be as large as A^pf < 30 - 80 22]. However, the spa- 
tial distribution of hadronic activity does not enter the 
dynamics of these simulations. It is conceivable that all 
hadronic activity, even if emerging from a large num- 
ber Ampi of partonic interactions, is located in few 'hot 
spots' in the transverse plane. Such a picture has been 
advocated previously in several contexts. For instance, 
models which view the proton as a collection of three 
black disks of diameter d result in a total p-p cross sec- 
tion, which matches experimental data for a surprisingly 
small diameter of d ~ 0.2 fm. Such simple constituent 
quark models account well for gross features of hadronic 
collisions, such as the ratios of hadronic cross sections 
(e.g. a-np/cjpp ~ 2/3) over a wide range of y/s. A more 
field theoretic motivation of a picture of hot spots may be 
based on the fact that a large value of Ampi requires high 
parton densities in the proton. In QCD evolution equa- 
tions such densities can arise from multiple branching of 
a few partonic components of high momentum fraction - 
and to the extent to which these branchings are collinear, 
partons within the proton wave function will be located 
indeed in few 'hot spots'. 

Irrespective of a specific dynamical picture, these con- 
siderations prompt us to consider scenarios for which the 
entire hadronic activity in a proton-proton collision is lo- 
calized in N s interactions regions of radius r$. Within 
each interaction region, density is distributed homoge- 
nous with Gaussian profile, and the regions are dis- 
tributed randomly according to the density profile of the 
proton (see left panel of Fig. [T] for N s = 3). To calcu- 
late the eccentricity for this class of models, we take into 
account that for finite A s the relevant eccentricity is the 
"participant eccentricity" [ll| 
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{y 2 }~{y} 2 , oxy = {xy}- 



{x} {y} and the event- by- event average {...} is taken over 
the distribution of interaction. This coincides with the 
definitions Eq. ([3]) in the limit of a homogeneous density. 
The resulting probability distribution for different sizes 
and different numbers of interaction regions is shown in 
Fig. [21 While the eccentricity vanishes indeed in the limit 




FIG. 2: Eccentricity distribution of central high multiplic- 
ity p-p collisions (central) for different number of interaction 
regions N s (top) of different size ro (bottom). 

of a homogeneous distribution (N s — > 00), we find that 
it is sizable for a wide parameter range. 

Elliptic flow develops on an event-by-event basis. How- 
ever, the experimental determination of U2 demands av- 
eraging over events and thus over distributions of geo- 
metric shapes. This averaging procedure is subtle; the 
relevant moment of the distribution P{e) is different for 
different methods of extracting V2 [23[ ■ We will focus on 
extracting v-i from 4-particle cumulants, «2{4}, for which 
e scales with [23J 
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Here, the average is performed over the -P(e) shown in 
Fig. [2] For A s = 3 our model yields the values of e{4} = 
0.69, 0.54, 0.40 for the different values of the spot size in 
a zero impact parameter collision. 

The precise transfer of spatial anisotropy e{4} to mo- 
mentum anisotropy t>2{4} requires a full dynamical cal- 
culation, but its magnitude can be estimated on general 
grounds [241 ]. For dilute systems, there is experimental 
evidence [8[ that the ratio v^jt depends linearly on the 
transverse matter density, 112/e oc (dN^jdy) / (S) [HI]. 

(S) 
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a 2 y denotes the mean trans- 



Here, (S) = 4tt^ 
verse area. For denser systems, the hydrodynamic limit 
is reached and the ratio becomes independent of the sys- 
tem size. In [2~ij ] the following interpolation formula was 
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FIG. 3: The flow signal 1)2 {4} as a function of multiplicity in 
most central p-p collisions, for models of N 3 — 3 interaction 
regions of radius ro = 0.1 (top 3 curves), ro = 0.2 (middle 
curves), and ro = 0.3 fm (bottom curves). Signals calculated 
iorX/Ko =0, 5.8 and 11.6 fm~ 2 are displayed by solid, dashed 
and dotted lines, respectively. 
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where (v 2 / ' e ^ hydro denotes the hydrodynamic limit, and 
K Q and A are constants which depend on the microscopic 
dynamics. The combination SX/ dN/dy can be under- 
stood as the Knudsen number of the system 24| . 

To estimate v 2 for p-p at the LHC, we assume that 
one can use Eq. |6]) with the input (v 2 /e) hvdro = 0.3 and 
X/Kq = 5.8 fm -2 extracted from the analysis of heavy 
ion collisions at RHIC [25j. This may be motivated by 
observing that many collective features depend mainly on 
event multiplicity and that the multiplicity of the p-p col- 
lisions considered here is comparable to that in peripheral 
lighter (Cu) ion collisions at RHIC, for which ([BJ applies. 
Under this assumption, we find for N s — 3 'hot spots' the 
large elliptic flow values v 2 {4} displayed in Fig. [3] Also 
for N s — 6 sufficiently small 'hot spots', the resulting 
v 2 {4} is found to be measurable (data not shown). To il- 
lustrate effects of theoretical uncertainties on deviations 
from equilibrium, we have varied in Fig. [3] the value 
of X/Kq between its hydrodynamic limit X/Kq and 
twice the fitted value quoted above. 

Not measuring v 2 in p-p collisions can have several 
reasons, such as a small spatial eccentricity of the in- 
teraction region, or viscous effects larger than expected 
from Eq. ([5]). However, combining the results shown in 
Fig. [3] and Table fl] we conclude that strong but phys- 
ically conceivable density fluctuations in the proton can 
lead to large values of the elliptic flow parameter v 2 in 
high multiplicity p-p collisions. Measuring elliptic flow in 
these collisions would provide a novel constraint on the 
number and distribution of multiple parton interactions 
within a p-p collision, which is a key input in the model- 



ing of the underlying event. At the same time, measur- 
ing flow in the smallest experimentally accessible system 
would show that collectivity develops on a sub-fcrmi time 
scale; since any fluid picture demands local equilibration, 
this measurement would provide a tight constraint on the 
equilibration and dissipation mechanisms of QCD. 

Note added. While preparing this letter, Ref. [26| ap- 
peared. By studying the eccentricity for a model of two 
flux tubes in a p-p collision, this work arrives at similar 
conclusions about the measurability of v 2 ■ 
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